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Abstract 

Laboratory wave basin measurements of the surge, heave and pitch of a floating plastic disk 
caused by regular incident waves are presented. The measurements are used to validate two the¬ 
oretical models: one based on slope-sliding theory and the other on combined potential-flow and 
thin-plate theories. 


1 Introduction 

As ocean surface waves progress deeper into the partially sea ice covered ocean, they encounter discrete, 


relatively thin chunks of ice (floes) of increasing horizontal dimensions (Squire & Moore, 1980). The 
range of dimensions depends on geographic location and season. However, the floes can be as small as 
a metre in the case of pancake ice, and up to hundreds of metres. The waves are attenuated by their 


interactions with the floes before they reach the quasi-continuous ice pack (Shen & Ackley, 1991), 


notwithstanding the large floes pushed into the outer fringes of the ice cover by random ice motions, 


which are subsequently broken up by the waves (Squire et al, 1995) 


The waves impact the ice cover. They breakup the ice into smaller floes (Prinsenberg & Peterson 


2011), which are more prone to melting and easily stirred up by winds, currents and waves. For 


example, waves herd the floes into groups (Wadhams 1983). Further, waves cause floes to collide with 


one another (Martin & Becker, 1987), which cause them to erode and produce rubble (McKenna Sz 


Croker, 1990). The collisions can turn into rafting events, which can lead to floes bonding and hence 


thickness growth (Dai et al., 

2004 

i. Waves also introduce warm water and overwash the floes, which 

accelerates melt (Wadhams et al . 

1979 

Massom & Stammerjohn 

2010 

)• 


Theoretical/numerical models have been developed to predict wave impacts on the ice cover. Shen 
& Ackley] (1991) used a one-dimensional model to study collisions between floes and herding. They 
used the slope-sliding model of Rumer et al. (1979) to calculate the horizontal motions of the floes 


induced by waves. The slope-sliding model is an extension of Morrison’s equation, which includes 
a force due to the slope of the wave field. The model is derived on the assumption that floes do 
not modify the wave field, i.e. the floe diameter is much less than the wavelength. It predicts the 
horizontal motion of a floe to be the sum of an oscillatory surge motion at the period of the incident 


wave, and a steady drift in the direction of the incident wave. Shen &; Zhong (2001) derived analytical 


solutions to the slope-sliding model in certain cases. Marchenko (1999) independently derived a similar 


identified an error in the derivation of Rumer et al. 


Kohout & Meylan (2008) and Williams et al. 


slope-sliding theory to Rumer et al. (1979). Grotmaack & Meylan (2006) related the two theories and 


(1979). 


(2013a b) modelled wave-induced breakup of a 


large group of ice floes. They applied breakup criteria that extended the earlier work of Langhorne 


et al. (2001). The kernel of both models is a model of a wave interacting with a solitary floe. The 
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wave-floe interaction model uses linear potential-flow theory to model water motions and thin-plate 
theory to model the floe. The linear potential-flow/thin-plate model is commonly used to study 
wave-floe interactions (see the review of Squire, 2007, for example). It assumes all motions are small 
perturbations from the equilibrium, i.e. the floe oscillates but does not drift. 

Kohout & Meylan (2008) and Williams et al. ( 201.3a|6 ) used two-dimensional models (one hor¬ 


izontal dimension and one depth dimension). Masson &; LeBlond (1989), Meylan et al. (|1997) and 


Bennetts et al. (2010) developed three-dimensional models of waves propagating through large groups 


of floes. They focussed on the attenuation of wave energy into the ice-covered ocean and did not model 


breakup or any other impact of the waves on the ice cover. Masson Sz LeBlond (1989) and Meylan et al. 


(1997) modelled the floes using the thin-disk models of Isaacson (1982) and Meylan & Squire (1996), 


respectively, noting the former is a rigid model and the latter is an elastic model. Bennetts et al. 


(2010) used a disk model and also square-plate model, using the finite-element approach of Meylan 
(2002), but did not find the different shapes significantly altered the attenuation rates predicted. 


Bennetts & Williams (2015) recently used laboratory wave basin experiments to validate the model 


of Meylan et al. (1997), and the two-dimensional model of Bennetts & Squire (2012), which was used 
by Williams et al. ( 2013a|6~ ). They used arrays of 40 to 80 identical wooden disks to model the ice 
cover, and measured the proportion of wave energy it transmitted for regular incident waves over a 
range of wave frequencies and, in some cases, for two different amplitudes. The quotient of thickness, 
r, over diameter, D, for the disks was t/D ~ 3.3 x 10~ 2 . The quotient of the incident wavelength, 
A, over the disk diameters was in the range X/D ~ 0.67 to 6.28. The incident steepness, represented 
by the product ka , where k = 2ir/X is the wavenumber and a is the incident amplitude, was in the 
range ka ~ 0.04 to 0.26. They showed the models predict the transmitted energy accurately for small 
incident amplitudes and low concentrations of the disks. They observed the models were inaccurate for 
the larger incident amplitudes when wave overwash of the disks was strong. (Overwash refers to the 
wave running over the top of the disks, due to their small freeboards.) Further, they provided evidence 
the models were inaccurate for high concentrations due to collisions between the disks, caused by out 
of phase surge motion of adjacent disks, and rafting, cause by out of phase heave and pitch motions. 
The potential-flow/thin-plate model does not include the highly nonlinear processes of overwash and 
collisions. 

Modelling collisions between disks requires an accurate model of the surge motion of a solitary disk. 
Heave and pitch motions must also be modelled to predict rafting. However, the potential-flow/thin- 
plate and slope-sliding model predictions of these oscillatory motions have not yet been thoroughly 
validated. In particular, the solitary disk experiences overwash for moderate incident amplitudes, 
which is not contained in either model. 


Bennetts & Williams (2015) presented measured surge, heave and pitch motions of a solitary 


wooden disk for a subset of the incident frequencies and amplitudes used for their multiple-disk tests. 
They compared the measurements to the predictions of the potential-flow/thin-plate model, and found 
the model is, in general, accurate. However, they also showed the model was least accurate for a test 
in which strong overwash occurred. In particular, the model overpredicted the translational motions, 
surge and heave, and underpredicted the rotational motion, pitch. 

Previously, Montiel et al. (20136 a) presented measurements of the oscillatory motions of a thin 
plastic disk in response to regular incident waves, as functions of the incident frequency. They used 
three thin disks, with t/D = 2.1 x 10~ 3 to 6.9 x 10 - ' 3 , and incident waves with lengths ranging from 
X/D ~ 0.63 to 3.14, and two small steepnesses ka ~ 0.03 and 0.06. They compared the measurements 
to predictions of the potential-flow/thin-plate model. However, they focussed on the flexural motion 
of the disk. They used a vertical rod through the centre of the disk to suppress surge, and a barrier 
around the edge of the disk to prevent overwash. 

Meylan et al. (2015) presented measurements of the surge, heave and pitch motions of a thin plastic 
disk, as functions of X/D. They used a disk with thickness over diameter quotient t/D ~ 3.8 x 10” 2 , 
and incident waves with lengths ranging from X/D ~ 0.9 to 12.3 and steepness ranging from ka ~ 0.01 
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to 0.3. They compared the surge measurements to predictions of the slope-sliding model. They 
showed the model predictions are accurate for incident wavelengths approximately greater than two 
floe diameters, for suitably chosen model parameters. However, they also used a barrier around the 
edge of the disk to prevent overwash. 


McGovern & Bai (2014) presented measurements of the heave and a composite surge and drift 
of model floes made of paraffin wax. They tested a variety of shapes, including square, rectangular 
and triangular shapes, but not disks. They modelled thick multiyear floes, and hence used relatively 
large thickness over characteristic length, D c , quotients, typically r/ D c = O(10 _1 ). They used regular 
incident waves with lengths in the range A /D c = 0.1 to 0.75, and steepness in the range ka ~ 0.03 
to 0.28. They studied heave and composite surge-drift as functions of A / D c and 2a/X, but did not 
compare these results to model predictions. They noted the occurrence of overwash for large incident 
amplitudes and high steepnesses, and suggested this may be the source of the reduced heave responses 
they found in this regime, which mirrors the finding of Bennetts <fc Williams (2015). 


McGovern & Bai (2014) also presented measurements of the drift of their model floes. They com¬ 
pared the measurements to the predictions of Stokes drift theory. They found the theory slightly 
underestimates the measurements. This finding is consistent with that of Huang et al. (2011). Lab¬ 


oratory experimental studies of the drift of floes have also been conducted by, for example, |Harms 
(1987) in a two-dimensional setting. 

In this study, surge, heave and pitch motions are extracted from an extended dataset to that used 
by Meylan et al. (2015). The extended dataset includes motions of a disk without an edge barrier, 
thus allowing inferences about the effect of overwash on the motions to be gained. The motions are 
presented as functions of A /D for three incident wave amplitudes. The motions for the disk without 
the barrier are shown to be consistently smaller than those with the barrier, especially for heave in 
the short-wavelength regime. Surge, which is the focus of this study due to its relevance to collisions, 
is also presented as a function of steepness, ka, for two incident wavelengths. The surge response is 
shown to decrease as steepness increases for the disk without a barrier and the shorter wavelength but 
is approximately constant otherwise. 

Further, the surge motions predicted by slope-sliding and potential-flow/thin-plate models are 
compared analytically and numerically in the long-wavelength regime. They are shown to have the 
same form. The surge, heave and pitch motion predictions of the potential-flow/thin-plate model are 
compared to the experimental measurements, as a function of A /D. In the short-wavelength regime, 
the model is shown to agree better with the data for the disk without a barrier. 


2 Experiments 

2.1 Method 

Laboratory experiments were performed at the Australian Maritime College, Launceston, Australia, 
using the Model Test Basin (mtb) facility. During the experiments, the wave induced motions of 
solitary floating disks were recorded. Figure [l] shows the plan view of the MTB and experimental 
set-up. 

The MTB is 35 m long and 12 m wide. It was filled with fresh water of density p =1000 kg/m 3 to 
a depth of h = 0.83 m. A piston-type wave maker bounds the MTB at its right-hand end. A sloping 
beach bounds the MTB at its left-hand end. 

Locations in the mtb are defined by the Cartesian coordinate system (x, y, z). The coordinate (x, y ) 
defines locations in the horizontal plane, parallel to the equilibrium water surface. The coordinate x 
points from the beach to the wave maker. The coordinate z defines the vertical location. It points 
upwards and its origin is set to coincide with the equilibrium water surface. 

Tests were conducted for a range of regular incident wave conditions. Target wave amplitudes from 
a = 2.5 mm to 50mm, and frequencies from / = 0.5 Hz to 2.0Hz, were tested. The corresponding 
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Figure 1: Schematic plan view of the mtb. Wave probes (Kl) record the wave profiles at set locations, 
and motion tracking cameras (■) record the motions of the disks. 
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/ = 1.5 Hz 
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Table 1: Summary of test matrices. 


wavelengths were approximately A = 0.4 m to 5 m. The wavelengths are calculated as A = 2it/h, where 
the wave number k is the positive real root of the dispersion relation 

uj 2 

fctanh kh = k where k = — (1) 

9 

is a frequency parameter, uj = 2irf is angular frequency, and g ~ 9.81 ms^ 2 is gravitational accelera¬ 
tion. 

Four wave probes were installed around the wave basin to record the generated wave properties. 
The measured incident frequencies closely matched the target values. The measured incident ampli¬ 
tudes were generally slightly smaller than the target amplitudes. The results presented in § |2.3| to §[4] 
use the target frequencies and the measured amplitudes. 

Tests were conducted for two matrices of incident wave amplitudes and frequencies. The first 
matrix contained more frequency entries than amplitude entries. The second matrix contained more 
amplitude entries. Table [T] summarises the tests conducted. Tests were not conducted for the largest 
amplitude and highest frequency combinations, to avoid wave breaking. 

Two thin plastic Nycel disks were installed in the mtb. Nycel is an expanded rigid foam PVC. The 
disks had radii R = 200 mm, thickness r =15 mm, density pd ~ 636kg m -3 , and hence equilibrium 
draft d ~ 9.55 mm and mass m = 1.2 kg. 

An edge barrier was attached to one disk, which is referred to as Disk B. The disk with no barrier 
is referred to as Disk NB. The left-hand panel of figure [2] shows a photo of Disk B. The edge barrier 
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Figure 2: Left-hand panel: photo of Disk B. Edge barrier, and rods and markers are visible. Right- 
hand panel: cross-sectional schematic of Disk B, which shows the mooring system. 



Figure 3: Schematic of Disk NB tests, including coordinate system and oscillatory motions. 


is a 50 mm high and 25 mm thick styrofoam ring. The barrier is used to prevent waves overwashing 
the surface of the disk, due to its small freeboard, and hence investigate whether the overwash affects 
disk motions. 

The locations of the centres of mass of the the disks at time t are denoted (x, y , z) = (X(t), Y(t), Z(t)). 
The two disks were initially positioned approximately halfway down the basin. They were placed 4 m 
apart to minimise their scattered waves interfering with one another. Similarly, they were each placed 
4 m from the MTB side walls, to minimise the effect of their scattered waves being reflected back to 
them. The origin of the Cartesian coordinate system in the horizontal frame is set to coincide with 
the geometric centre of Disk B in the horizontal plane, and, thus, (A'(0), Y (0)) = (0, 0) for Disk B and 
(A(0),T(0)) = (0,4) for Disk NB. The disks were anchored to the floor via loose elastic tethers to 
prevent them from drifting too far down the tank and to assist in resetting the initial positions after 
each test. 

Four light-weight tracking balls were attached to each disk via aluminium rods. Figure [2] shows the 
markers and rods on Disk B. The position of each marker was captured by the Qualisys non-contact 
motion tracking system. The system consists of eight pairs of infrared cameras and receivers installed 
along the perimeter of the wave basin. Qualisys records the locations of each marker in real-time at 
200 frames per second. It uses this information to calculate the translational and rotational motions of 
the disks. The markers were placed at different heights to minimise the chance of them overlapping in 
the camera image. Data was collected from t = 0, when the wave maker was activated, up to t = 60 s, 
when the wave maker was turned off. 

Motions in the six rigid degrees of freedom of the disks were recorded. The symmetric motions, 
with respect to the x-axis, were dominant due to the plane incident waves. The translational motion in 
the x-direction is a combination of drift, restoration due to the mooring system, and surge, which is an 
oscillatory motion at the frequency of the incident waves. The translational motion in the ^-direction 
is oscillatory heave. The rotational motion in the x^-plane is oscillatory pitch. Figure [3] illustrates the 
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Figure 4: Left-hand panel: example motion in ^-direction. Steady-state interval is denoted by dashed 
vertical lines. Right-hand panel: close-up of smoothed signal (solid curve). Local maxima and minima 
are denoted by circles. Mean maxima and minima are denoted by dot-dashed lines. 


oscillatory motions, which are analysed in this investigation. 

2.2 Data processing 

Figure [4] shows an example time series provided by Qualisys for the translational motion of Disk B 
in the ^-direction. The test shown used an incident wave with frequency / = 1.25 Hz and measured 
amplitude a = 8.5 mm. 

Motions are considered in the steady-state interval only, which begins after the initial transient 
phase of the motions pass, and ends when waves reflected by the beach begin to interfere with the 
motions. The left-hand panel of figure [4] denotes the intervals with vertical dotted lines. The interval 
is calculated using the phase velocity of the wave, c = /A. 

The MATLAB smooth function with the lowess method is used to eliminate noise from the raw 
signal and clearly identify the local maxima and minima. The smooth function applies a local regression 
method using a weighted linear least squares algorithm. The degree of smoothing is controlled by the 
smoothing parameter — a factor of 0.01 is specified here. The right-hand panel of figure [4] shows the 
smoothed signal for a 5 s interval at the beginning of the steady-state interval. 

The maxima and minima are identified, as shown by the circles in the right-hand-panel of figure 
[|| The heave amplitude, ap, is calculated as half the difference between the average peak and trough 
values. The same procedure is used to calculate the pitch amplitude, ap, from the time series of the 
translational motion. 

The left-hand panel of figure [5] shows the translational motion in the x'-direction, for the same test 
considered in figure |4j Surge and drift are visible in the motion, as is the long-period restoration due 
to the mooring system. 

The smooth function is used with a smoothing factor of 0.9 to extract the drift from the Qualysis 
data. The surge motion is the difference between the full data and the drift. The right-hand panel 
of figure [5] illustrates the decomposition for the small time interval. The surge amplitude, as, is, 
subsequently, calculated using the same procedure as for heave and pitch. 
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Figure 5: Left-hand panel: as in left-hand panel of figure [4] but for motion x-direction. Right-hand 
panel: close-up of smoothed motion (solid curve), and decomposed into surge (dot-dashed) and drift 
(dashed). Local maxima and minima of surge motion are denoted by circles. 
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Figure 6: Raos for Disk NB, as functions of nondimensional incident wavelength. Data are grouped 
according to incident wave amplitude: a =10mm (□), 20nun (o) and 40mm (A). 
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Figure 7: Example of changes to the wave held in the long- and short-wavelength regime (left-hand 
and right-hand panels, respectively). Scattering is evident in the short-wavelength regime. 
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Figure 8: Raos for Disk B (x) and Disk NB (o). 


2.3 Response amplitude operators 

Figure [6] shows response amplitude operators (raos) of Disk NB for surge, heave and pitch, as functions 
of the incident wavelength nondimensionalised with respect to the disk diameter, D = 2 R. The results 
are for the first test matrix. The surge RAO is as/{acoth kh}, where acothkh is the horizontal extent 
of the trajectory of a fluid particle at the free surface. The heave RAO is ap/a. The pitch RAO is 
ap/ka , where ka is the incident wave steepness. In the sub-panels, data are grouped according to 
incident wave amplitude (using different symbols and colours). Error bars are omitted here for clarity. 
They are presented in figures [13] to [15] 

The RAOs are independent of the incident wave amplitude. The mean range of the corresponding 
RAOs for the different amplitudes is approximately 0.067. Moreover, no behavioural trend with varying 
amplitude is evident. The results suggest a linear relationship exists between the amplitudes of 
oscillatory motion and the incident wave amplitude. 

The RAOs are approximately unit value for X/D greater than 3. The disk does not affect the 
incident waves in this long-wavelength regime. The left-hand photo in figure [7] shows an example of 
this behaviour. The disk simply follows the path of a fluid particle at the free surface. 

The RAOs decrease as the incident wavelength decreases. The disks scatter the incident waves in 
the short-wavelength regime. The right-hand photo of figure [7] shows an example of this. Scattering 
results in less energy being transferred into the oscillatory motions, and thus reduces the RAOs. 

Figure [8] compares the RAOs of the two disks for the tests in the first matrix. Tests with different 
incident amplitudes are not distinguished from one another. Disk NB generally has slightly smaller 
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Figure 9: Deviation of surge RAO from the mean, as a function of wave steepness, for incident frequency 
/ = 1.25 (left-hand panel) and 1.5 (right). Results for Disk B (x) and Disk NB (o). Error bars indicate 
deviation of the maxima/minima. 


RAOs than Disk B. The difference is greatest for heave at 14.8% and least for surge at 7.8%. The 
percentages quoted represent the mean of the difference in mean RAO for the two disks, for an incident 
wavelength and RAO combination. Differences are most pronounced in the short-wavelength regime, 
X/D < 3. The largest difference is 0.35, which occurs for heave at X/D = 0.9. 

Figure [9] shows the surge RAO, as a function of the incident wave steepness, ka, calculated from 
tests in the second matrix. Surge motion is presented here, as it is the focus of this study. Cognate 
behaviours occur for heave and pitch motions (not shown). 

The results are grouped according to the incident wavelength used (different panels) and which 
disk the results refer to (different symbols and colours). The RAOs are presented in terms of their 
relative deviation from the mean, defined by 

. Surge RAO — (Surge RAO) 

A RAO = — r 

(Surge RAO) 


( 2 ) 


where angled brackets (•) denote the mean with respect to both disks and different steepnesses. The 
error bars denote RAO values for the maximum and minimum surge amplitudes in the tests. 

For the longest incident wavelength tested, X/D = 2.50, the RAO values generally deviate by 
no more than 10% from the mean. The deviations here do not display a consistent trend and are, 
therefore, attributed to measurement errors. For the shortest wavelength tested, X/D = 1.73, the 
RAOs for Disk B are insensitive to steepness. In comparison, the RAOs for Disk NB decrease as 
steepness increases. This is attributed to the wave overwashing the Disk NB in this regime, and hence 
suppressing surge. This finding is consistent with that of Bennetts & Williams (2015). 


3 Theoretical Models 

3.1 Slope-sliding model 

The slope-sliding model assumes the floating body, here the disk, moves along the wave profile due to 
gravity. Its movement is resisted by drag between the body and water. The equation of motion in the 
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x-direction is thus 


x dV 

m(l + c m ) — = -mg 


dg 

dx 


+ pc d l\V\V. 


(3) 


x=X 


This nonlinear ordinary differential equation is solved for the horizontal velocity of the disk, V(t) = 
dX{t)/dt. 

The term on the left-hand side of equation ([3]) is the inertial force of the disk. The coefficient c m 
is the added mass of the disk, i.e. its increased resistance to motion due to contact with water. The 
first term on the right-hand side is the sliding force due to gravity. The quantity r/(x, t ) is the wave 
profile. The second term is the drag force. The quantity V = V w (t) — V (t), where V w (t) is the velocity 
of water particle on the free surface at x = X. The coefficient Cd is the drag coefficient and l = ttR 2 
is the wetted surface area of the disk. 

The model further assumes the disk diameter is small in relation to the incident wavelength. As 
discussed in 


2.3, in this regime the incident wave profile is not modified by the floe, i.e. 


r](x,t ) = a sin(A;x — cut). 


(4) 


The phase of the incident wave does not affect the RAOs extracted from the model, and it is therefore 
normalised to the origin. 

Using the wave profile Q and the notation Q = kX — cut, equation ([3]) reduces to the autonomous 
dynamical system 

dV —mgka cos Q + pcji IVIV 

- —, (5a) 


dt 


where V(V, Q ) = \cuasm.Q — V, and 


m( 1 + Cm) 


d Q , T/ 

dt 


(5b) 


Figure [TO] shows the phase plane for an example problem, with zero drag and added mass, Cd = c m = 0, 
and an incident wave of amplitude a = 50mm and length A = 6m. 

System ([5]) always has a limit cycle. Figure 10 illustrates the limit cycle with the solid (red) curve 
that intersects the origin. The limit cycle represents steady surge (the amplitude of the cycle) and drift 
(the difference in period of the cycle and the incident wave period) of the disk. This is the relevant 
solution for the present investigation. 

The example problem also possesses two fixed points at (Q,V) = (7r/2,c o/k) and (37r/2, co/k). 
These are the so-called surfing solutions. They can only be reached for a large initial velocity, V(0). 
The right-hand fixed point is a centre, i.e. the trajectories surrounding it are closed. The left-hand fixed 
point is a saddle node, which has lower and upper homoclinic connections. For all initial conditions 
below its lower homoclinic connection in the phase plane, the orbits tend to the limit cycle. 

For small non-zero values of the drag coefficient, the fixed points are shifted and the right-hand 
fixed point becomes a stable spiral. However, the system retains the structure that all orbits below 
the lower homoclinic connection tend to the limit cycle. For large values of the drag coefficient, the 
fixed points disappear in a saddle-node bifurcation. All orbits then tend to the limit cycle. 

In practice, the limit cycle is obtained by solving system ([5]) as the second-order ordinary differential 
equation 

d 2 X 


m(l + Cm) 2 + mgka cos(kX — cot ) — pcdW\V\V = 0, 


( 6 ) 


with zero initial displacement and velocity, i.e. X = dX/dt = 0 for t = 0. The solution to ([6]) is 
obtained numerically using the MATLAB package ode45 solver, which is based on a fourth and fifth- 
order Runge-Kutta method. The steady solution, i.e. the limit cycle, is decomposed into surge and 
drift components, as for the experimental data. 

The values of the added mass and drag coefficients are obtained via comparison with experimental 
data. Meylan et al. (2015) compared the surge RAO obtained from the slope-sliding model to the 
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Figure 10: Phase diagram for horizontal disk velocity, V, and phase, Q, for example problem with 
u = 2.68, c m = Cd = 0, a = 0.05 m and A = 6 m. Vectors show the evolution of V and Q. Here two 
equilibrium points exist at (1.57, 2.56) and (4.71, 2.56). The lower red curve represents the limit cycle. 


experimental surge data for Disk B. They found the parameters c m = 0.1 and Cd = 0 provided a good 
fit to the data for A /D > 2, approximately. For shorter incident wavelengths, the slope-sliding model 
does not capture the surge for any choice of added mass and drag coefficients. 


3.2 Linear potential-flow/thin-plate model 

The linear potential-flow/thin-plate model assumes the disk oscillates about its equilibrium position. 
The oscillations are time harmonic at the frequency of the incident wave, and are symmetric with 
respect to the x-axis, i.e. in-line with the incident wave. The location of the centre of mass of the floe 
in the x-direction is denoted 

X(t) = R e{A s e~ iut }, (7) 

where As is the complex-valued surge amplitude, such that its magnitude is the surge amplitude, as, 
and its argument is the phase of surge motion. The vertical location of its lower surface is denoted 
Z(x,t) = — d + Re{?n(x)e^ lwt }. The displacement function w is decomposed as 


w(x) = A h + A P x + ^2 £ nw n {x )• 


( 8 ) 


71=1 


Here Ap and Ap are complex-valued amplitudes of heave and pitch motions, respectively, such that 
ap = \Ah\ and ap = \Ap\ . (Here, pitch is defined with respect to rotations around the geometric 
centre of the disk’s lower surface, as opposed to the centre of mass. This has a negligible effect for the 
thin disk considered.) 

The functions w n are the (symmetric) flexural modes of vibration and are the corresponding 


amplitudes. Itao &; Crandall (1979) and Meylan & Squire (1996) provide expressions for the flexural 


-5 


mm, 


modes. For the incident wave and floe parameters studied here, the amplitudes £n < 5 x 10 
i.e. the disk responds to the incident waves in its rigid modes only. 

Disk motions are forced by differential pressures from the surrounding water due to hydrodynamics. 
Air pressure is assumed to be constant, po say. Water pressure is modelled using the linearised version 
of Bernoulli’s equation 

< 9 $ 


P = Po - P 


dt 


pgz. 


(9) 
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Here $ is the velocity potential of the water, i.e. the water velocity field is the spatial gradient of <f>. 
Following potential-flow theory, it has been assumed the water is inviscid, homogeneous, incompressible 
and in irrotational motion. 

Water motions are also assumed to be time harmonic, and the velocity potential is expressed 


$(x,y,z,t) = Re\-^-4>(x,y,z)e la; *) . 

^ 1 CO ) 


( 10 ) 


The reduced (time independent) velocity potential, 4>, satisfies Laplace’s equation throughout the 
water domain, i.e. 


V 2 0+^-^ = O where V = (d/dx, d/dy), 
oz A 

and a no-normal-flow condition on the floor of the basin, i.e. 

d<P , 

—— = 0 on 2 = —n. 
oz 


( 11 ) 


( 12 ) 


The water is assumed to extend to infinity in all horizontal directions, i.e. reflections from the basin 
boundaries are not considered. In the far-held the velocity potential is composed of the incident-wave 
potential 

ae lkx cosh k(z + h ) 


4>i{x,y,z) = 


(13) 


cosh kh 

and a geometrically decaying scattered-wave potential, which satisfies the Sommerfeld Radiation con¬ 
dition 

y/r (d r — ik) (</> — 4>i) —>• 0 as r —> oo, (14) 

where r = y/{x 2 + y 2 } is the radial coordinate. 

The amplitudes of the waves and disk motions are assumed to be sufficiently small that linear theory 
is applicable. Conditions on moving boundaries are therefore approximated by linearised conditions 
applied on the corresponding equilibrium boundaries. The linearised free-surface condition is 


dcj) 

dz 


= K(fi on z = 0, 


(15) 


which holds at all horizontal points (x, y) ^ H, where = {x, y : x 2 + y 2 = r 2 < R 2 } is the projection 
of the hoe in equilibrium onto the xy-plane. 

Following Kirchhoff-Love thin-plate theory (Timoshenko & Woinowsky-Krieger 1959), the lin¬ 
earised equation of motion in the z-direction is 


(1 — k d)w + FV^w = (f> on z = —d, 


(16) 


and for (x,y) E H. Here F = Eh?/{Ylpg(\ — iz 2 )} is a scaled hexural rigidity of the disk, where 
E = 530 MPa is the Young’s modulus, as measured by a three-point bending test on a strip of the 
Nycel. The quantity v = 0.3 is chosen as a typical value of Poisson’s ratio for Nycel. Note, the model 
is insensitive to the elastic parameters, E and u , as the disk responds rigidly. 

The motion of the disk and the velocity potential are also coupled via the linearised kinematic 
conditions 

d0 


— = kcos( 0)X on (x,y) E Sfi, —d < z < 0, 
or 


where 6 = tan (y/x) is the azimuthal coordinate, and 

dcj) 


dz 


= kw on (x, y) E H, z = —d. 


(17a) 


(17b) 


12 






Lastly, the vertical displacements of the floe satisfy the linearised free-edge conditions 


V^-(l-zz) -—+ 


1 dw 1 d 2 


w 


r dr 


r 2 QQ2 


d 


Id fid 2 


+ (!-«/)-- - 


dr 


w 


r dr \r d9 2 


0 and 

(18a) 

= o, 

(18b) 


both for (x, y) E SO. 

Equations ([nJi-Q define a boundary value problem for the velocity potential, </>. The vertical 
displacement of the floe, w, is obtained as part of the solution. At this juncture, the surge amplitude, 
Ag, is a parameter of the problem. 

The velocity potential and vertical displacement function are decomposed as 


<f> = 4> + Agcj) and w = w + Agw. 


(19) 


The functions <f> and w are the solutions of equations (11) to (18) with the floe artificially restrained 


in surge, Ag = 0. The functions (j) and w are the solutions of equations (11) to (18) with no incident 


wave forcing, a = 0, but with the disk forced to oscillate with unit amplitude in surge, Ag = 1. These 
boundary value problems are solved using an eigenfunction matching method, which is described by 


Linton & Mclver (2001) for general hydrodynamic problems. 


The surge amplitude, Ag. for the full problem is found via the linearised equation of motion of the 
disk in the x-direction, which is 

- (w 2 m + Cm)Ag = f. (20) 

The surge motion is forced by the pressure field around its edge created by the incident wave, i.e. 

/ = P9 [ [ Rcos{9)[(t>\ {Xty)£ sn d9 dz. (21) 


' —d J — it 


The complex quantity C m is defined by 


/ 0 rn 

/ Rcos(9)[j)}( Xty ) & sn d-9 dz. 
-d J —it 


( 22 ) 


The imaginary and real components of C m are known as the added mass and damping terms, re¬ 


(19). 


spectively (Mei, 1983). The full solutions, <and w, are obtained by substituting Ag into equation 


3.3 Long-wavelength regime 


The slope-sliding model is valid in the long-wavelength regime. The left-hand panel of figure 11 
shows the drift velocity predicted by the slope-sliding model tends to zero as the incident wavelength 
becomes large. The figure shows this holds for zero and non-zero drag and added mass. The slope- 
sliding model, therefore, predicts the horizontal motion of the disk is dominated by oscillatory surge 
in the large-wavelength regime. 


For zero drag, which gives the best fit to the experimental data (Meylan et al. 2015), the slope- 
sliding model is 

dy 


,dV 

(1 + Cm) ^ — 9 


dx 


= —gakcosfkX — uit) 


(23) 


J x=X 


The small drift is translated into the condition \kX\ <C 1. The slope-sliding model is, therefore, 
manipulated into the form 

(1 + c m) < -^r = ~ gakcosujt + O(kX), (24) 


13 























Figure 11: Left-hand panel: drift velocity versus non-dimensional wavelength, as derived from the 
slope-sliding model. Here c m = Cd = 0 (solid black curve), c m = 0 and Cd = 0.5 (black dashed), 
c m = 0.1 and Cd = 0 (red dot-dashed), c m = 0.1 and c m = 0.5 (dotted red curve). Right-hand panel: 
motion in the x-direction, calculated from the slope-sliding model, using the parameters A /D = 15, 
and c m = Cd = 0. The surge component of the full numerical solution is represented by the blue 
dashed curve. Analytical long-wavelength solution is denoted by the solid black curve. 


which is integrated to give 


V' ( t)=- g f 8 iD (‘f) + 0 ( fcX). 


w(l + C m ) 

Thus, to leading order in kX, the displacement of the disk is 


X(t) = 


gak cos (ut) a coth(kh ) cos(o;t) 
oj 2 (1 T c m ) 1 T c m 


(25) 


(26) 


where the final expression is derived using the dispersion relation. The slope-sliding model, therefore, 
predicts the surge amplitude as = a coth(kh) /(1 + c m ) in the long-wavelength regime. The right-hand 
panel of figure [XT] shows the leading-order displacement (26) accurately approximates the surge motion 
of the full solution for X/D = 15 and c m = Cd = 0. 

The long-wavelength solution for the potential-flow model can also be solved analytically by con¬ 
sidering the horizontal equation of motion in equations (20) to (22). In long waves, the velocity 


potential is approximately equal to the incident potential (j)j, as with equation (13), due to the wave 
field remaining unaltered by the floe. 

In (13), the term e lkx = e lkRcosd can be expressed as 1 + ikRcosO + O(kR) 2 through Taylor 


expansion. Since kR < 1 in long waves, O(kR) 2 may be disregarded. Substituting the remaining 


terms into equation (21) gives 


/ = PS [ [ aR cos 9(1 + ikRcosO )—-— 'f — —- d9 dz « iagk nR 2 dp (27) 

J-dJ-TT cosh kh 


as k tends to zero (i.e. in long waves). Note that irR 2 dp = m as required by Archimedes’ principle. 
Equating this with the left-hand side of (20) gives 


A S = 


—i agkm 
U) 2 m + Crr. 


(28) 
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X/D X/D 

Figure 12: Added mass and damping coefficients (left and right-hand panels, respectively), calculated 
from the potential-flow model, as a function of non-dimensional wavelength. 


hence the surge amplitude is 


as = 14s I = 


a coth kh 

11 + C m I 


(29) 


after applying the dispersion relation. In the above, C m = C p +iC a where C p is the damping coefficient 
and C a is the added mass coefficient nondimensionalised with respect to u 2 m. These two coefficients 


can be calculated using equation (J22|) . 
Figure 


12 


shows the value of the real and imaginary parts of C m for increasing wavelengths, 
calculated using dimensions of the disks described in §|2.1| In the long-wavelength limit, both the added 


mass and damping coefficients tend to zero. At the largest wavelengths considered, i.e. X/D = 12.3, 
C a = C p = 0(1(T 3 ). Based on this, the assumption is that C m ~ 0 in long waves, and therefore the 
surge amplitude will tend to as = a coth kh. 


4 RAOs: potential-flow/thin-plate model and data comparison 

The left-hand panels of figures [l3] to [l5] show the RAOs predicted of the potential-flow/thin-plate 
model for surge, heave and pitch, respectively, as functions of non dimensional incident wavelength. 
Predictions are shown for the disk model outlined in § 3.2 and a simplified two-dimensional model. 


The model predictions are overlaid on the RAOs calculated from the experimental data, as given in 
figure [8} Error bars are included for the experimental results here. As in figure [9j the upper and lower 
limits are calculated using maximum and minimum amplitudes for each test. The right-hand panels 
show the corresponding differences between the data and the disk model. 

Visually, the potential-flow/thin-plate models predict the RAOs measured during the experiments 


with pleasing accuracy. The right-hand panels of figures 13 to 15 show the differences between the 
disk model and the experimental data. Here, experimental data are divided into 14 groups, with each 
group corresponding to a particular wavelength. The means for each group are plotted (separately for 
Disks B and NB) in terms of its difference from the disk model. Error bars on the right-hand panels 
represent the overall maxima and minima for each group. Differences between the means and the disk 
model are generally very small (less than 30%) across the range of wavelengths considered. For Disk 
B, the average difference is greatest for heave at 6.8%, and least for surge at 1.2%. For Disk NB, the 


15 












Figure 13: Left-hand panel: comparison between potential-flow models (disk: solid curve, 2-D: dashed) 
and experimental data for surge as per figure [8j Error bars show the range of error for each data point. 
Right-hand panel: experimental data for each disk are grouped according to wavelengths, and the 
means for each group are plotted in terms of its difference (d RA o) from the disk disk model. Error bars 
represent the overall maxima and minima of each group. The comparison with the 2-D potential-flow 
model (dashed) is also shown. 
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Figure 14: As per figure [l3| but for heave. 
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Figure 15: As per figure [I3| but for pitch. 


average difference is greatest for heave and pitch at -5.6%, and least for surge at -3.9%. The general 
trend 

Differences between the two-dimensional and disk models are evident only for X/D < 10, approx¬ 
imately. The differences are generally very small across the range of wavelengths considered. The 
differences are greatest for surge at 4%, and least for pitch at 1.1%. However, the two-dimensional 
model predicts a phase change in surge at X/D ~ 6, which is not predicted by the disk model. This 
implies the models will be dissimilar for shorter wavelengths. 

Figure [IB] shows the deviation between the disk model and the experimental data from the second 
test matrix. The deviation is presented in a similar fashion to figure [9j but with the model predic¬ 
tions replacing the mean, The model generally underpredicts the measured surge RAOs for the longer 
wavelength, X/D = 2.50. There are two exceptions, both for Disk NB. The model also underpredicts 
the surge motion of Disk B for the shorter wavelength, X/D = 1.78. The model underpredicts the 
surge of Disk B for the lowest steepness only. The decrease in surge RAOs for Disk NB with increasing 
steepness results in the model overpredicting the surge motion. 

5 Summary and conclusions 

Wave basin experiments were conducted to determine the hydrodynamic responses of an ice floe in 
regular waves, over a range of wave conditions. Wave amplitudes between 2.5 and 50mm, and wave 
frequencies between 0.5 and 2.0 Hz were tested. Scale models of typical small-medium sized floes were 
used. The floes, which were constructed from Nycel, measured 400 mm in diameter and 15 mm in 
thickness. An edge barrier was used to determine the effect of overwashing. Results for the single floe 
experiments showed that: (a) surge, heave and pitch motions are, in general, linear with respect to the 
incident wave amplitude; (b) there was a reduction in all three motions in the short-wavelength regime 
as a result of scattering; (c) overwashing occurs in the short-wavelength regime and is dependent on 
the incident wave amplitude; (d) overwashing appears to decrease the response of all three motions 
— the greatest reduction occurred in heave, while surge was least affected. 

Two theoretical models were presented — a slope-sliding model and a linear potential-flow/thin- 
plate model. Comparisons with experimental data showed that the potential-flow/thin-plate model 
provided good agreement across the full range of wavelengths tested. The slope-sliding model was 
only accurate in the long wavelength regime. 
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Figure 16: As per figure [9j but here, deviation is calculated with respect to the potential-flow model 
(horizontal solid line). 

A simplified 2-D potential-flow model was also proposed, by disregarding the lateral dimension. 
Although there were minimal differences in the responses predicted by the simplified and full models, 
the latter had a slightly better agreement with the experimental data. The two-dimensional model, 
however, was still considered to be reasonably accurate. 
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